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Abstract 

We investigate the exact bright and dark sohtary wave solutions of an effective one dimensional (ID) Bose-Einstein 
condensate (BEC) by assuming that the interaction energy is much less than the kinetic energy in the transverse 
direction. In particular, following the earlier works in the literature Perez-Garci'a et al. [50], Serkin et al. [51], 
Giirses [52] and Kundu [53], we point out that the effective ID equation resulting from the Gross-Pitaevskii (GP) 
equation can be transformed into the standard soliton (bright/dark) possessing, completely integrable ID nonlinear 
Schrodinger (NLS) equation by effecting a change of variables of the coordinates and the wave function. We consider 
both confining and expulsive harmonic trap potentials separately and treat the atomic scattering length, gain/loss 
term and trap frequency as the experimental control parameters by modulating them as a function of time. In the 
case when the trap frequency is kept constant, we show the existence of different kinds of soliton solutions, such as the 
periodic oscillating solitons, collapse and revival of condensate, snake- like solitons, stable solitons, soliton growth and 
decay and formation of two-soliton like bound state, as the atomic scattering length and gain/loss term are varied. 
However when the trap frequency is also modulated, we show the phenomena of collapse and revival of two-soliton 
like bound state formation of the condensate for double modulated periodic potential and bright and dark solitons 
for step-wise modulated potentials. 
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1. Introduction 



Experimental realization of trapped Bose-Einstein condensates (BEC) in alkali metal atoms has triggered 
immense interest in understanding the various properties of ultra eold matter [1,2]. These include four wave 
mixing, formation of vortices, bright solitons, dark solitons, gap solitons, interference patterns and domain 
walls in binary BECs [3-14] . An interesting dynamical feature in the context of Bose-Einstein condensation 
is the formation of matter wave solitons [4-12,15-19]. For example, dark solitons [7-11] and gap solitons [12] 
were experimentally observed in BECs with repulsive interactions, whereas bright solitons were demonstrated 
in systems with attractive interactions [4, 5] . The recent experiments at Heidelberg and Hamburg universities 
have shown the formation of dark solitons, their oscillations and interaction in single component BECs of 

Rb atoms with confining harmonic potential [20-22]. These experimental results have been augmented 
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by extensive theoretical work including predictions on bright and dark matter wave solitons [15-17], lattice 
solitons [18], and soliton trains [19]. 

The dynamics of a BEC at absolute zero temperature is usually described by the mean-field Gross- 
Pitaevskii (GP) equation [23], which is a generalized form of the ubiquitous NLS equation [24-26], for the 
wave function of the condensate. The NLS equation and its generalizations are widely used in the nonlinear 
optics literature [24-26] to describe the evolution of coherent light in a nonlinear Kerr medium, envelope 
dynamics of quasi-monochromatic plane wave propagation in a weakly dispersive medium, high intensity 
pulse propagation in optical fibers and so on. 

As is well known the NLS equation in one space and one time dimension is a completely integrable soliton 
system [24-26]: It possesses bright solitons in the focusing case, while dark solitons exist in the defocusing 
case. On the other hand, the GP equation which is a three space and one time dimensional nonlinear 
evolution equation, is not integrable in general, but can admit exact solutions only in very special cases. In 
general, the GP equation has to be analyzed numerically or by approximation methods. For the cigar-shaped 
trap potential when the transverse linear oscillator length is much smaller than the longitudinal length, the 
3D GP equation can be reduced to the effective ID GP equation by assuming the kinetic energy of the 
longitudinal excitations and the two-body interaction energy of the atoms are of the same order, while both 
are much less than the kinetic energy of the transverse excitations [27,28]. To obtain bright and dark soliton 
solutions, its quasi-lD version has to be analyzed essentially numerically. 

Interestingly, some of the soliton solutions in the presence of regular harmonic confinement, with appro- 
priately tailored gain or loss, have been shown to exhibit collapse and revival phenomena, with an increase 
in amplitude. For example, Hulet and his co-workers observed the growth of a condensate of trapped ''Li 
atoms with attractive interaction and subsequent collapse [29], and Miesner et al. measured the rate of 
growth of a ^^Na condensate fed by a thermal cloud [30]. Such a growth of a BEC from thermal vapor 
was also experimentally realized for atoms of *^Rb [31]. Several methods have been introduced theoretically 
to account for the growth of BEC [32-35]. The GP equation with the complex linear gain or loss term, 
often referred to as GP gain equation, was used to model this type of BEC dynamics by Drummond and 
Kheruntsyan [35]. It has been shown that, as the condensate grows, the center of mass oscillates in the trap. 
The GP gain equation is of special importance in the field of atom laser, which is a device that produces an 
intense coherent beam of atoms by stimulated process [36,37]. In recent years, several attempts have been 
made to understand the dynamics of BEC by studying the GP gain equation [38-40] . 

Further, there has been an increased interest during recent times in studying the properties of BECs 
with time varying control parameters. In particular, investigations have been made over a decade or so to 
understand the properties of condensates under 

(i) the temporal variation of atomic scattering length which can be achieved through Feshbach reso- 
nance [41-44], 

(ii) inclusion of appropriate time-dependent gain or loss terms which can be phenomenologically incorpo- 
rated to account for the interaction of atomic cloud or thermal cloud [40] , 

(iii) periodic modulation of trap frequencies [45] and so on. 

Earlier studies on BECs have shown that a suitable periodic modulation of scattering length can stabi- 
lize the condensate against collapse [46-48]. However, there are certain occasions where a sign alternating 
nonlinearity and increase in the frequency of oscillations result in the acceleration of collapse [49] . In this 
connection, a systematic study on the condensates under the influence of time-dependent parameters which 
can be achievable through experiments is of considerable importance. 

Most of the theoretical studies on the matter wave solitons in BECs have been carried out using either 
numerical or certain approximation (perturbation) methods [15-17]. For example, Ruprecht et al [15] have 
shown the existence of bright soliton solutions in weakly interacting ID and 2D BECs by solving the 
corresponding GP equations numerically. In an another context, Pcrez-Garcfa et al [16] have shown the 
existence of bright soliton solutions to the GP equation with negative scattering length in highly asymmetric 
traps by using the multiple-scale approximation method, while Busch and Anglin [17] have shown the motion 
of dark solitons in trapped effective ID BEC by using a multiple time scale boundary layer theory. On the 
other hand, there are only a few exact analytical methods available in the literature [40,50-53] to deal with the 
GP equation. In this connection it is worth mentioning that Serkin et al [51] have constructed exact soliton 
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solutions by relating the effective ID GP equation to the compatibility condition of a Lax pair which is in turn 
related to that of the Lax pair of the ID integrable NLS equation using an inverse scattering transform (1ST) 
method while Atre ct al [40] transformed the GP equation to the unforced Duffing equation and thereby 
derived soliton solutions from the elliptic function solutions of the Duffing equation. Also, Perez-Garcia 
et al. [50] obtain a similarity transformation connecting the NLS equation with time varying coefficients and 
the autonomous cubic NLS equation. In addition, there exists several other interesting works in the literature 
to construct exact soliton solutions of the ID GP equation: Specifically, for certain time-dependent atomic 
scattering lengths (nonlinearity) with (i) time independent homogeneous, linear and quadratic potentials 
without gain/loss term [54], (ii) time independent expulsive parabolic potential without gain/loss term [55], 
(iii) time independent expulsive parabolic potential and arbitrary constant gain/loss term [56] , (iv) tanh type 
modulated harmonic trap frequency without gain/loss term [57], explicit growing/decaying type solutions 
have been constructed. In the optical soliton context, in Ref. [58] exponentially growing solutions have been 
constructed for either constant atomic scattering length or zero potential case with a specific tanh type 
dissipation term. However, a more detailed analysis for bright and dark soliton solutions in quasi-lD BEG 
with time-dependent atomic scattering length, harmonic trap frequency and gain/loss term is of considerable 
present interest from an experimental point of view. 

In this paper we bring out a large class of bright and dark soliton solutions by simply mapping the effective 
ID GP equation onto the standard NLS equation, following the earlier works in the literature [50-53]. In 
particular we show that the effective ID equation resulting from the Gross-Pitaevskii (GP) equation in 
confining as well as expulsive potentials with time varying parameters, namely the scattering length, gain/loss 
term and trap frequency, can be transformed into the standard soliton (bright/dark) possessing, completely 
integrable ID NLS equation by effecting a change of variables of the coordinates and the wave function 
under certain condition, which is related to the scattering length, gain/loss term and the trap frequency. We 
consider both confining and expulsive harmonic trap potentials separately and treat the atomic scattering 
length, gain/loss term and trap frequency as experimental control parameters by modulating them as a 
function of time. We have deduced different kinds of bright / dark soliton solutions for different forms of the 
trap frequency by suitably tailoring the atomic scattering length and atom gain/loss term such as 

(i) periodic oscillating solitons, collapse and revival of condensate, and snake-like solitons for time inde- 
pendent confining harmonic potential, 

(ii) stable solitons, soliton growth and decay, and formation of two-soliton bound state for time independent 
expulsive harmonic potential, 

(iii) Oscillatory solitons, phenomena of collapse and revival of two-soliton like bound state formation of the 
condensate for double modulated periodic potential and 

(iv) stable solitons, soliton growth and decay for step-wise modulated potential. 

We also make comparisons of our results with the experiments of Khaykovich et al [5] and Strecker et al [4] 
obtained for ^Li atoms and for ^^Rb atoms [20-22]. Further, we also compare our results with the theoretical 
results of Atre et al [40] and Serkin et al [51]. 

It may be further noted that in experiments the matter wave solitons [4, 5, 20-22] have been observed 
within few micrometer (/im) lengths and times up to the order of few ms. In the present paper we choose the 
soliton parameters such that the length (/im) and time (ms) scales are well within the experimental regime. 
Further, we have confined ourselves to the issue of stabilizing solitons by choosing the forms of modulated 
atomic scattering length and atom gain/loss term so as to balance each other. We have also considered 
here only the two-body interactions and neglected the three-body recombination loss which will introduce 
a quintic term in the GP equation which does not in general reduce to exactly integrable standard NLS 
equation. Such a case may be studied by the numerical and approximation methods. Another interesting 
aspect of the temporal modulations is to look into the problem of parametric amplification, which however 
falls outside the scope of the present analysis and should be dealt with separately. It may also be noted 
that in the context of optical systems with varying coefficients, only picosecond optical solitons have been 
observed [59-65]. 

This paper is organized as follows. In Sec. 2, we discuss the method of transforming cylindrically symmetric 
GP equation with time varying parameters such as the atomic scattering length, gain/loss and with both 
confining and expulsive harmonic potentials into an effective ID GP equation with time varying parameters. 
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Then, wc map the ID GP equation to the standard NLS equation by making use of suitable transformations, 
subject to a constraint on the forms of scattering length, gain/loss term and the trap frequency. Next, in 
Sec. 3 we deduce the generalized bright and dark matter wave soliton solutions of the GP equation from 
the soliton solutions of standard NLS equation. We show the existence of various forms of bright and dark 
soliton solutions, their collapse and revival, and decay and growth for suitable choices of the condensate 
parameters with time independent confining and expulsive harmonic potentials in Sees. 4 and 5, respectively. 
This is followed by a demonstration of various forms of bright and dark matter wave soliton solutions for 
temporally modulated potentials in Sec. 6. Finally, in Sec. 7, we present a summary of the results and 
conclusions. In Appendix A, we have shown the validity criteria of the ID GP equation for various forms 
of control parameters studied in this paper. In Appendix B, we give a brief description on the solution 
procedure of solving the restrictive condition for which the effective ID GP equation can be mapped on to 
the standard NLS equation. We also tabulate the possible solvable choices of the modulated trap frequency 
forms from the Ricatti equation which specifies the restrictive condition in Appendix C and provide a table 
of physical meaning of symbols used in the text in Appendix D. 



2. Gross-Pitaevskii equation w^ith complex gain or loss term 



In the present study, we consider the GP equation with the inclusion of a complex gain or loss term of 
the form 



ih 



9^'(r,t) 
dt 



2m 2 



(1) 



where U — ^^^as{t), as{t) is the time varying s-wave scattering length, T4rap(r, i) = Vo{x,y) + Vi{z,t), 
Vo{x,y) = Y'^J-i^'^ + y"^) is ^ cylindrical time independent harmonic trap, Vi{z,t) = ^^uj1[t)z'^ is a time- 
dependent harmonic trap along the z direction which can be either confining or expulsive, ct = ±1 depending 
on whether the interaction is attractive (cr = +1) or repulsive {a = —1), and T{t) is the gain/loss term which 
is phenomenologically incorporated to account for the interaction of atomic cloud or thermal cloud. The plus 
sign of the T{t) term leads to the mechanism of loading external atoms (thermal clouds) into the BEG by 
optical pumping while the minus sign of T{t) describes a BEG that is continuously depleted (loss) of atoms. 
If the condensate is fed by a surrounding thermal cloud, then the condensate undergoes a periodic growth 
and collapse. The above equation has gained much importance in recent times to understand the growth 
and decay of the condensate due to time-dependent parameters [38-40]. However, one has to note that the 
parameters are varied sufficiently slowly so as to satisfy the adiabaticity criteria so that the wave-packet 
continues to evolve in the same quantum eigenstate as a function of time. A possible condition for the 
validity of the adiabaticity criteria [66,67] is that the characteristic time scale, say T, of the variation of the 
modulated parameters is slow enough when compared with the instantaneous nonlinear eigenvalue E-a\{t) of 
the nonlinear equation at time t, such that E^\{t)T / {2t:) ^ 1. 

We also note that for the GP equation (1) the normalization condition in the presence of the gain/loss 
term r{t) becomes 



f l^-pd^r = TVexp / Tit)dt 

J —OO \_J 



(2) 



where N is number of the condensate atoms. 

There are typical characteristic length scales to specify the condensate density. The confining frequencies 
of the harmonic trap potential (Vtrap), that is lo±_ and u}z{t), set the length scale for spacial size which are the 
transversive oscillator length a± = ^^'^ longitudinal length az{t) = ( ntJ'(t) )"^^^- addition, 

the effective mean- field nonlinearity introduces the length scale, namely the healing length ( ~ 8TTnas{t), 
which is the distance over which the kinetic energy and the interaction energy balance [28,68,69]. Here 
n = N / {a\az{t)) is the mean density of the condensate. 

The presence of the nonlinear term in the above GP equation (1) complicates its solution procedure. 
One should note that it is very difficult to obtain closed form solution of Eq. (1) as such. In most cases the 
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above equation can be studied either by approximation methods or by suitable numerical methods. However, 
under certain circumstances, Eq. (1) can be reduced to a simpler form by appropriate transformation whose 
solutions are known. One way here is to reduce or relate the above Eq. (1) to the standard NLS equation in 
(1+1) dimensions so that the soliton solutions can be readily obtained. In the following we shall describe a 
method of reducing Eq. (1) to the NLS equation and obtaining the consequent spatiotemporal patterns. 

First let us reduce the GP equation (1) to an effective ID equation which deals with the cigar-shaped trap 
potential when the transverse linear oscillator length a± satisfies the following criteria: a± is much smaller 
then the longitudinal length az{t) [28,68,69], that is 



< 1, a_L ~ C- (3) 



The above criteria essentially means that only for an axially open trap {a± ^ <iz{t)) one obtains a ID GP 
equation. In this case we assume that the condensate wave function can be taken as 

*(r, i) = - cxp (-lij^t - \^ \ $( A^^^i). (4) 

Here the width and length of the condensate are of the order of a±_ and a^, respectively, and hence n = 
N /{a\az) is the mean density of the condensate. Note that the kinetic energy of the longitudinal excitations 
fi^ /2ma1{t)] and the the two-body interaction energy of atoms (/i^as/m)n] are of the same order, 
while both are much less than the kinetic energy of the transverse excitations h^/2ma\) [27,28]. From 
this, we get following criteria for the validity of the reduction of the 3D GP equation to the effective ID GP 
equation in the absence of the gain/loss term T{t) [27,28], 

^ « 1. (5) 

az[t) 

However, in the presence of the gain term in Eq. (1), and in view of the normalization condition (2), the 
above criterion may be modified by replacing N with A^exp [Jr(t)(ii] as 

2a,(t)jVexp [jT{t)dt\ ^ ^ 
az{t) 

Note that due to the above condition, for a given N, the functional form of the scattering length as(t), 
gain/loss term V{t) and the longitudinal(axial) length az{t) are to be chosen appropriately as a function 
of time so that the above condition can be valid for sufficiently long times, which we will demonstrate in 
Appendix A for various forms as(t), T(t) and az{t) studied in this paper. 

Substituting the above form of the wave function (4) in Eq. (1), one obtains the effective ID GP equation 
with modulational gain/loss term, atomic scattering length and trap frequency as 

i^, = -1$,, - aR{t)mH + + z^<i>, (7) 

where f^^ = 7 = = = and as is the Bohr radius. Note that in eq (7), 

the variable z actually represents z/a±, due to the transformation (4). Similarly t stands for uj±t. Using the 
relation (4), we can rewrite the above condition as 



1^12,3 2NaB 
|$|^d'^r = exp 



j{t)dt 



(8) 



In addition, one can also transform away the gain (or loss) term ^{t) in Eq. (7) by applying the following 
transformation [70] 



$(z,<) = exp 



Q{z,t), (9) 



5 



so that Eq. (7) can be rewritten in terms of the new variable Q{z,t) as 

iQt^-^Qzz~'jR{t)\Q\^Q+^^z^Q, (10) 

where R{t) is a time-dependent parameter and can be written as a function of R{t) and ^(t) in the following 
way 



R{t) = exp 

2.1. Transformation to standard NLS equation 



-f{t)dt 



R{t). (11) 



It has been shown by Giirses [52,53] recently that the above Eq. (10) can be mapped onto the standard 
NLS equation under the following, transformation, 

Q = Ag(Z,r), (12) 

where the new independent variables T and Z are chosen as functions of the old independent variables t 
and z as 

T=G{t), Z = F{z,t), (13) 

while A = A(2;,i) is a function of t and z.It may also noted that similar mapping has been done by Serkin 
et al. [51] and Perez-Garcia et al. [50] using the lens transformation. Applying the above transformation 
(12), so that 

^-gA+F,^ ^-fA ^-f ^ + f-^ (G.-^F-A) (14) 

dt~ 'dT^ 'dz' dz~^dz' dz^~ ''dz^^dz^' ^ '~ dt'^~ dz' ^ ' 

one can reduce Eq. (10) to the standard NLS equation of the form 

«9t + 9zz + crlgj^? = 0, (15) 
subject to the conditions that the functions A, F, G, and R should satisfy the following set of equations, 

«At + ^A,, - ^z^A = 0, (16a) 

iKFt + i (2A,K + AF,,) = 0, (16b) 

Gt = ^^ (16c) 

The above NLS equation (15) is widely used in the nonlinear optics to describe the evolution of coherent 
light in dispersive medium, high intensity pulse propagation in optical fibers and so on, which is completely 
integrable equation and admits soliton, elliptic and trigonometric function solutions [24-26] . 

In order to solve the unknown functions A, F, G in the above equations (16) we assume the polar form 

A = r(z,i)exp[i6'(z,t)], (17) 

One can immediately check from the relations (16c) that r is a function of t only, r = r{t), since G and R are 
functions of t only. Then from Eqs. (16) one can easily deduce the transformation function A given by (17) 
and the transformations G{t) and F{z,t) through the following relations (Details are given in Appendix B): 

=== 2rlR, (18a) 

9 = + 2brlRz - 2b^rf^ [ R^dt, (18b) 

2R J 



F{z,t) = 2rQRz - 46rg J RMt, (18c) 
G(t) = 2rl [ R^dt, (18d) 
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Here b, tq are arbitrary constants, and R and il^ should be related by the following condition (see Ap- 
pendix B, Eq. (B.12) for made details) 



0, 



(19) 



which is a Riccati type equation with dependent variable and independent variable t. Note that here 



R(t) = R{t)exp 



exp 



7(i)di 



is a quantity representing the modulated scattering 



length and gain/loss term and f2(t) — LOz{t)/uj± represents the time-dependent trap frequency. Eq. (19) 
implies a specific constraint between these quantities and for a suitable choice of them satisfying the equa- 
tion (19), the effective ID GP equation (7) and so Eq. (10) can be mapped onto the standard NLS equa- 
tion (15). It may be noted that a transformation similar to the above transformation has also been reported 
by Serkin et al. [51] by relating Eq. (10) to the compatibility condition of the related Lax pair and Perez- 
Garci'a et al. [50] using the lens transformation. However, we will make use of the relations (12)-(19) as they 
are more straightforward. Also we note that the transformation given in Eq. (18) is more general than the 
one given in Refs. [50,51] because of the additional term proportional to the arbitrary constant b in Eqs. (18 
b) and (18 c). 



2.2. Solution to the Riccati equation (19) 



One can consider three possible approaches to the solutions to the present Riccati equation constraint (19): 

(i) Given 51^(t), obtain possible explicit forms of R{t) by solving (19), 

(ii) Given R{t), deduce the form of the trap frequency Q,{t) as can be written down from (19) straightfor- 
wardly 

(iii) Given ri^(i), obtain numerical solution for R{t) when the Riccati equation cannot be explicitly solved 
for R{t) 

Firstly, there are numerous forms of r2^(i) for which explicit solutions for R{t) can be obtained [see for 
example Ref. [71]] by solving the Riccati equation (19). Most important of them are given in the form of 
Table C.l in Appendix C. Next, for a given form of R{t), the corresponding form of the trap frequency Vl{t) 
can be straightaway calculated from (19) as 



\ 




(20) 



Finally, for those forms oiVL^it) for which no exact solutions can be found by solving the Riccati equation (19), 

numerical solutions can always be found for from which R can be calculated. One can make use of all 
the above three approaches to obtain the matter wave solutions. 



3. Bright and dark soliton solutions 

The NLS equation (15) has been known for quite some time for its novel soliton solutions. It admits 
Jacobian elliptic function solutions, bright and dark solitary wave solutions and trigonometric function 
solutions [24,72], depending on the sign of a. For instance, soliton solutions of the standard NLS equation (15) 
were given by Zakharov and Shabat by solving the Cauchy initial value problem through the 1ST method 
for (7 = -|-1 (focusing or attractive case) [73] and for cr ~ —1 [74] (defocusing or repulsive case). The 
corresponding soliton solutions are called bright (cr = +1) and dark (a = —1) solitons, respectively. In 
the following, we will make use of the one soliton solutions only and higher order soliton solutions will be 
discussed separately. We obtain the one soliton solutions of the effective ID GP equation (7), by making use 
of the one soliton solutions of NLS equation (15) through the relations (9), (12) and (18). 
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3.1. Attractive (focusing) type nonlinearity (a = +1) 



The bright sohton solution of (15) for the case of a 
unimportant phase constants) 



-1 takes the following form [74] (leaving aside the 



{Z, T) =aexp 



2 



2a' - c' 



-T 



sech 



(21) 



Here the real parameters a and c correspond to the amplitude and velocity, respectively, of the bright 
soliton envelope and plus sign in the superscript represents the a = +1 case. Correspondingly for the GP 
equation(7), using the relations (9), (12) and (21), the bright soliton solution can be written as 



$+(2, t) = exp 



lit) 



dt 



'{z, t) = exp 



lit) 



dt 



A(z,t)g+(Z,T). 



(22) 



Now substituting the forms of A(z,t), Z and T from Eqs. (13), (17) and (18) in the above Eq. (22), we get 
a generalized expression for the bright soliton as 



$+(z,t) ^A+{t) sech [^+{z,t)] exp [ir]+{z,t)] , 
Here the time-dependent amplitude of the generalized soliton (23a) is given by 

A+{t) =aroJ2R{t)eI = aro^/2i^(^)e/ 



(23a) 



(23b) 

which depends on the form of the atomic scattering length Rit) and gain/loss term jit), and the generalized 
wave variable which specifics the width and velocity of the bright soliton as a function of time and 
longitunumberdinal distance, is given by 



^+{z,t) =V2aroRz- V2arl{c + 2brn) J R^dt = V2araR 



z-ro(c + 2&ro)i / R^'dt 
R 



(24a) 



so that the width is proportional to l/\jR(t) and the velocity is a function of time. Also the generalized 
phase of the bright soliton is given by the expression 



7J+{z,t) 



-R 



^z'' + roic + 2bro)Rz-r^ 



2R 



(c + 26ro)2 



R'dt. 



(24b) 



Note that the amplitude of the generalized bright soliton ^^(t) given by Eq.(23b) can vary as a function 
of time depending upon the forms of the scattering length R{t) and gain/loss term 7(i). Consequently one 
can have growth or decay or constancy (of the amplitude) of the bright soliton for appropriate choices of Rit) 
and jit). In particular for the specific choice R{t) = exp[— 2 j j{t)dt], the amplitude = constant and the 
soliton amplitude remains constant in time, which has not been noted in the earlier literature [40,51,55-58]. 

However note that the corresponding width (cx; l/\fR = / '^^^^'^^ / R{t)) as given (24a) need not be 
constant which can increase or decrease as the case may be. In the following sections 4-6, we will discuss 
difi^erent kinds of interesting generalized soliton solutions for different forms of unmodulated/modulated 
harmonic potentials by suitably tailoring the atomic scattering length and gain/loss term. 

Finally the normalization condition (8) can be rewritten in terms of the soliton parameters as 



/ |$+(z,t)|2dz = 2%/2aroCxp ( 

J —00 J 



j{t)dt 



2NaB 



■ exp 



j{i)di 



which allows one to fix the constant rg in terms of the experimental BEC parameters as 

Nub 



ro 



V2a 



(25) 



(26) 
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3.2. Repulsive (defocusing) type nonlinearity [a = — 1) 



On the other hand, when cr = — 1, one may write the dark one sohton solution in the form [73] 



(Z, T) = ^ [c - 2i/3 tanh /3(Z - cT)] exp 
v2 



(27) 



where the parameters /3, c correspond to the depth and velocity of the dark soliton, respectively. Here the 
minus sign in the superscript means the cr = — 1 case. 

Using Eqs.(9), (12) and (27), one can write down the dark soliton solution of (7) as 



$ (2, t) = exp 



7(t) 



dt 



'(z, t) = exp 



lit) 



dt 



Aiz,t)q-{Z,T) 



Now substitute A(z,t), Z and T from Eqs. (13), (17) and (18) in the above Eq. (28), and we get 

7(i). " 



$ {z,t) = ^roV i?exp 



-dt 



{c — 2i/3tanh (z, t)] } exp [177 {z,t)] 



(28) 



(29a) 



Here the generalized wave variable ^ which specifies the width and velocity of the dark soliton as a function 
of time and longitudinal distance is given by the expression 



C iz,t) = 2(3roRz - 2/3r^ (c + 26ro) J R^dt = 2/3roi? 
and the generalized wave phase of the dark soliton is expressed as 



-ro(c + 2fero)i / R'dt 
R 



T]-(z,t) ^ iz^ + 2br?.Rz^ 

2R 



rl (c^ + 2b^ 



+ 413^) / R^dt 



Also the depth of the generalized dark soliton is of the form 

lit) 



= { 2/3roV^(Oexp 



-dt 



I = |2/3ro/RW 



exp 



-/{t)dt 



(29b) 



(29c) 



(30) 



which depends on R{t) and jit). As in the case of the bright soliton solution we can get stable dark 
soliton solutions by choosing the forms of R{t) and ^{t) to balance each other so as to achieve a constant 
amplitude/depth. In this case also, the width depends on i?(i) and j(t) and can increase/decrease or remain 
constant with time. Again the normalization condition (8) of the dark soliton solution for the GP equation 
leads us to the condition 

2NaB 



4Ric' + A(3')~\^-{z,t)\' 



dz = 4/3ro exp 



j(t)dt 



■ exp 



"fit)dt 



so that the parameter tq can be fixed in this case as 



NaB 
2l3a±' 



(31) 



(32) 



3.3. Choice of condensate parameters 



Now making use of the forms of R for given 51^ (i), as presented in Table C.l in Appendix C or through 
a detailed numerical analysis, into the above bright or dark soliton solutions (23a) or (29a), respectively, 
one can analyze the nature of the BEG matter wave properties. In Table G.l, we list out a chosen set of 
27 different choices of ft{t) and the corresponding forms of R{t) for which the Riccati type equation (19) 
exhibits closed form solutions. In the following sections, we shall analyse the nature of the bright and dark 
BEG solitons, their growth, decay, collapse and revival for a select set of physically interesting choices of 
n{t) and R{t) from Table G.l. 

Let us first focus our attention on the case of time- independent harmonic potential for which 0^(t) in (7) 
is a constant (51^ (t) — In this case, we shall identify the different soliton solutions for both confining. 
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fig > 0, as well as expulsive, Hq < 0, potentials for various forms of scattering length R{t) and gain/loss 
term j(t). It is of relevance here to point out that in recent experiments on the formation of bright matter 
wave solitons in an ultracold ^Li gas, both the confining as well as expulsive potentials have been used [4,5]. 
In addition, such a formation of matter wave bright soliton has been achieved by temporal modulation of 
atomic scattering length due to Feshbach resonance [4,5]. Further, as pointed out earlier, we note that the 
formation of dark solitons, their oscillations and interaction in single component BECs of Rb atoms with 
confining harmonic potential have been demonstrated in recent experiments at Heidelberg and Hamburg 
universities [20-22]. Motivated by these experiments, we fix the parameters such as the trap frequencies 
(w_L and ujz), number of atoms (TV) and scattering length in the transverse direction, a±, to confirm to the 
experimental values. In particular we choose the parameters as ui± ~ 2tt X 710Hz, iV = 2 X 10**, and a±_ = 1.4 
/im. The trap frequency in the axial direction is chosen as Wz = 27r x 70Hz for the case of the confining 
potential, while we choose lOz = ^iri x 70Hz for the expulsive potential. Correspondingly fig takes the value 

UJ-, 2 

JIq = — - « 0.01 so that rio is taken as 0.1 in our analysis in the following. Accordingly, in the entire 

manuscript the unit of time is taken as milliseconds and the length is expressed in micrometers. Also it is 
not difficult to provide a simple estimation of the validity of the reduction of 3D GP equation to the ID GP 
equation given by Eq. (5) in terms of dimensional units. For example, when one considers a BEG of Rb 
atoms with N « 10^, for Oz w 4.4 /im Eq. (5) gives the estimate as{t) < 0.22nm, which is experimentally 
realizable. Note that, in the presence of the gain/loss term the criteria for the validity of ID approximation 

becomes (see Eq. (6)) exp [J j{t)dt\ ^ 1. In the present analysis, we suitably tune R{t), ft{t) 

and "f{t) to satisfy the above criterion, which may be still realizable (which will be shown in the Appendix A 
for various forms). 

In the next section we shall discuss the formation of both bright and dark matter wave solitons in time 
independent confining harmonic potential for different choices of temporal modulation of the atomic scat- 
tering length and gain/loss term. This is followed by a demonstration of the solitons in time independent 
expulsive harmonic potential case in Sec. 5. Finally, In Sec. 6, we discuss the formation of both bright and 
dark matter wave solitons in time-dependent confining potential such as double modulated periodic potential 
and step-wise modulated potential. 



4. Time independent confining harmonic potential 



We consider the time independent confining potential with the ratio of trap frequencies obeying the 
relation r2^(<) = lip, where fto is a constant. Using 51^ (i) = JIq in Eq. (19), we obtain a particular solution 
for R as, [see Table G.l] 

R = sec{not + 5). (33) 

Since R = R{t) exp(/ j{t)dt), one can identify different types of soliton solutions of Eq. (7) for appropriate 
forms of R and 7. Accordingly, we will choose these forms and investigate the nature of matter waves. 



4.1. Case 1:^ = 0, bright soliton with periodic oscillations in amplitude 



When the gain term is absent in Eq. (7) so that for 7 = the atomic scattering length can be written, using 
Eqs. (11) and (33), as i? = sec(J7oi + S)- Accordingly, substituting the form (33) into (23a) the magnitude 
of the bright soliton solution of (7) is given by 



[$+(^,01 



NoE 



■ ^1 sec{^ot + S)\ > sech 



Nob 
a± 



scc{flot + S) \ z 



tan(f2ot + 5) 



V2q 



(34) 



It may appear that the presence of sec(r2ot + 5) in the the amplitude of the above bright soliton will lead to 
singularities at finite time. However, at the same time when sec(i7o^ + (5) — > 00 the envelope term [ sech(...)] 
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also goes to zero sufficiently fast so as to keep the overall amplitude finite. Note that the width of the 
soliton is proportional to cos(J7oi + 5)\ and hence it is varying periodically with time. Fig. 1(a) depicts 
the spatiotemporal pattern of condensate wave function in a confining oscillator potential for a typical choice 
of parameters. It is clear from Fig. 1(a) that the bright matter wave soliton wave for R = sec{Q.Qt + 5) and 
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Fig. 1. (a) Bright soliton for the case of Q'^(t) = SIq, 7 = 0, and R{t) = scc{Qot + S) and (b) the corresponding contour plot. 
The parameters are chosen as a = 1, Qo = 0.1, 5 = 0, b = 0.1, and c = 0.2 in Eq. (34). 

7 = shows periodic oscillation in amplitude due to the presence of the periodic terms sec(Ooi + S) and 
tan(57oi + ^) in the right hand side of Eq. (34). Fig. 1(b) shows the periodic oscillation of the width. 
In this case the corresponding dark soliton solution can be written as 

Nag 



\^-{z,t)\ 



2l3aj_ 



\/\se(iiht + 6)\ 



X < + 4/3^ tanh^ 



Nob 

aj_ 



scc{flot + S) \ z 



tan(riot + S) 

^0 



(35) 



However, in the absence of a suitable gain/loss term (7 = 0), the above expression becomes singular at finite 
values of t and so is not physically interesting. 

4.2. Case 2: 7 = ait — 02 sin(Kt), The collapse and revival of bright/dark solitons 



A physically interesting case can be realized by choosing the gain term as ^(t) = ait — 02 sin(Kt), (ai, a2 
and K are parameters), so that 



R{t) — scc{rtot + S) cxp 



ait^ a2 cos(Ki) 



(36) 



The above form of 'y{t) was suggested recently by Atre, Panigrahi and Agarwal [40]. In this case, the bright 
soliton solution of (7) can be written as 



Na 



|$+(z,i)| =^^V|sec(17ot + (5)|.exp 



a± 
X sech 



ait^ 02 cos{Kt) 



Nat 



sec(riot + 5) \ z 



\/2aa 



no 



(37) 



The amplitude of the soliton is -^/j sec(flQt + S)\. exp ^ — h °^ '^2^ ^-'^'^ ^he width is proportional 

to cos{flot + 5)\, indicating collapse and revival of the condensate. Fig. 2 illustrates the collapse and 
revival of the condensate wave function with an overall increase in the amplitude. The collapse and revival 
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Fig. 2. (a) Bright matter wave soliton for r2^(t) = Qq, 7 = ait — a2 sin(Kt), R(t) = scc{Cot + (5) with a = y^^^^, Qi = -001, 
a2 = 8.25, 6 = 0, c = 0, 5 = 0, f2o = 0.05, k = 1.5 in Eq. (37) and (b) the corresponding contour plot. Note the periodic 
collapse and revival of the wave with increase in amplitude over time. 

of atomic condensate and amplification [see Fig. 2] through periodic exchange of atoms with the background 
is obviously due to the partially sinusoidal nature of the gain function, 7 ~ a\l — a-i sin(Kt). 

We also note here that, Strecker ct al. have experimentally observed the formation of the bright soliton of 
^Li atom in quasi-lD optical trap, by magnetically tuning the interactions in a stable BEC from repulsive 
to attractive [4] . Theoretically, Atre, Panigrahi and Agarwal have shown that by suitably tailoring the gain 
profile in the same range of the experimental parameters, the matter wave soliton exhibits dramatic collapse 
and revival of condensate with an increase in the amplitude [40] . 

The corresponding expression for dark soliton can be written as 



NaB 
2/3a_L 



\J\ sec(r2oi + (5)1 exp 



02 COs(Kt) 



2k 



4/3^ tanh^ 



Nob 



sec(fioi + f^) ) ^ ^ 



2/3ai' 



tan(riot + S) 
1^ 



(38) 



The above expression, as in the case of the bright soliton (37), clearly shows that a periodic collapse and 
revival of the matter wave dark soliton can also be realized for the choice of R(t) as given by (36). The dark 





t (in ms) 



Fig. 3. (a) Dark matter wave soliton for r2^(i) 



a\t — a2 sin(ftf), R{t) = sec(f2of + 5) with /3 = -y^^ 



ai 



.001, 



a2 = 8.25, f) = 0, c = 0, 5 = 0, f2o = 0.05, k = 1.5 in Eq. (38) and (b) the corresponding contour plot, again portraying a 
collapse and revival with increase in amplitude over time. 



soliton (38) is depicted in Fig. 3. 
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4.3. Case 3: 7 — — fig tan(iloi + f^); snake-like bright/dark solitons 



Finally wc consider the case with another form of gain term as 7 = — f2otan(f2oi + S) so that R = 
sec^lflot + S] which can be obtained from Eqs. (11) and (33). In this case, the bright soliton solution of (7) 
can be written as 



|$+(z,OI =^ sech 
a± 



NOE 

a± 



sec(fioi + S) ] z — 



(39) 



One may note that in this case the amplitude of the soliton (A^ 



NaE 



) remains constant in time. However, 



the corresponding width cos(r2o^ + i5)|) is varying periodically in time, leading to a snake-like effect. 

Figs. 4 shows the bright solitary wave (39) exhibiting a snake-like effect. It is interesting to point out that a 
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Fig. 4. (a) Bright soliton for fl'^(t) = Slg and 7 = —Qq tan{Qoi + = sec[Oo* + (b) The corresponding contour plot. 

The other parameters are a = 1, 6 = 0.1, c = 0.5, <5 = 0, and f!o = 0.1 in Eq. (39). 

similar snake-like effect has been demonstrated by Serkin and Hasegawa in the context of dispersion managed 
optical soliton with sinusoidally varying dispersion coefficient, nonlinearity, and gain or absorption term [38]. 




The dark soliton solution, in this case, is given by 
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|$-(z,i)l 



NaB 
'2l3a± 



4/3^ tanh^ 



sec(rioi + 5) ] z — 



(40) 



As pointed out earlier in Sec. 4.3 for the bright sohton, the above dark soUton (with constant amphtude and 
periodically varying width) also exhibits a snake-like effect. Fig. 5 depicts the propagation of snake-like dark 
solitary wave. One may note that a similar type of snake-like dark matter wave solitons and their collisions 
have been observed in the recent experiments of '^^Rb condensates [20-22] with confining potential. In this 
context, it is quite interesting to have an exact solution of the form (40) for the snake-like dark soliton 
solutions for confining potential with time varying scattering length and atom gain/loss. 



5. Expulsive potential (Jl'^(t) — —^q is constant) 



In this section, we consider the case of time-independent expulsive potential il'^{t) = — f^o- ^'^^ ^l^is 
choice, substituting Q'^{t) ~ in the integrability condition (19), one may identify three independent 
particular solutions of physical interest for R [71], namely (i) R = exp(— fio^), (h) R = exp(f2oO: ^"^^ (iii) 
R = sech(51oi + ^) [see Table C.l]. However, the stabilization of bright and dark matter wave solitons 

is solely dependent on the constancy of the product term \l R(i) exp[J '^^^dt\ = ^ R(t) cxp[J ^{t)dt]. The 
corresponding restriction on the parameters is called safe range parameters. When the parameters are chosen 
outside the safe range, the bright and dark solitons may become unstable, that is, one may witness either 
growth or decay, as studied for (i) 7 = and R{t) — e^"* in Ref. [55] and (ii) 7 = constant and R{t) = 
sech(iloO 6"*^* ill I^sf. [56]. In fact, stable solitons can be formed only if we choose the form of the scattering 
length R{t) and the gain/loss term j{t) as R{t) = exp[— 2 J j{t)dt], so that the soliton amplitude remains 
constant in time. In the following, we discuss various kinds of soliton solutions of Eq. (7) for different forms 
of -f{t) and R{t) both within as well as outside the safe range parameter values in all the above three cases 
of R. 



5.1. Case 1: R = cxp(rioi) 



For R = exp(f2o<) case, we study the solution of Eq. (7) for both the negative (loss) and positive (gain) 
values of 7. In this case, we find stable solitons for the safe range value 7 = — fio and growth of matter wave 
soliton for the out of safe range value 7 = JIq, where J^o > 0. 



5.1.1. 7 = — Oq, stable bright/dark soliton: 

First let us consider the inclusion of depletion/loss term in (7), that is, 7 = — f2o < 0. In this case R = 
exp[2r2oi]: obtained from Eqs. (11) and (33), implying that the atomic scattering length is exponentially 
growing with time. Accordingly the bright soliton solution of (7) is given by 



|$+(z,i)| 



NOE 



secli 



NoE 



exp(not) ) z — 



^aa\ 



exp(2not) 



217o 



(41) 



The solitary wave solution occurs due to the balance between the negative gain term and the exponential 
growth of the atomic scattering length. It may be noted that, here, the overall amplitude remains constant 
(A"'" = '^^^). Fig. 6(a) shows the bright solitary wave (41) for 7 = — ilo and R = exp[2f2oi]. However, the 



width of the wave (ex 



~ exp(— irJoO) decreases exponentially with time as seen from Fig. 6(b). 



For the safe range parameter 7 
expression 



-fio the corresponding dark soliton solution of (7) is given by the 



|$-(z,0| = 



Nan 
2^ax 



+ 4/32 tanh^ 



( NaB 



exp(f2oi) ) ^ — 



N^al 



2l5'^a\ 2Pa] 



exp(2r2o<) 
2^0 



(42) 



One may again note that the amplitude remains constant (A = ) in the above expression. Figs. 7 
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Fig. 6. (a)Bright soliton for Q'^{t) = — t^Q and 7 = ~flo, R(t) = e^°'. (b) Corresponding contour plot. Tlic other parameters 
are a = 1, b = 0.5, c = 1.0,r2o = 0.1 in Eq. (41). 
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Fig. 7. (a) Dark soliton for Sl^(t) = —fig and 7 = —Ho, R(t) = c^"*. (b) Corresponding contour plot. The other parameters 
are /3 = 1/^2, 6 = 0.5, c = 1.0, Qq = 0.1 in Eq. (42). 

illustrates the dynamics of the dark solitary wave. Once again, as in the case of bright soliton, the width of 
the wave is decreasing exponentially with time [see Fig. 7(b)], while the amplitude remains constant. 

5.1.2. 7 = f2o; growth of bright/dark soliton: 

For positive gain term 7 = > case, one finds that R is independent of time, which indicates that the 
atomic scattering length is constant. In this case, the bright soliton solution of (7) can be written as 



|$+(^,i)| =:^^ sech 



Nof 



cxp{flot) 1 z — 



N'^a% \ exp(2rJo<) 



\/2aa^ 



2Qn 



exp [ftot] 



(43) 



One may note that, in the above expression, the amplitude of bright soliton grows exponentially (A^ = 
-^^^ exp(r2ot)) due to the presence of the gain term (parameters in the out of safe region), while the width 
(oc exp(— iript)) decreases exponentially with time. Fig. 8(a) depicts the exponential growth of bright solitary 
wave for 7 = JIq and R = constant. Here the width of the wave decreases with respect to time exponentially 
as seen in Fig. 8(b). 

In this case the dark soliton solution of (7) can be written as 
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Fig. 8. (a) Bright soliton for n^(t) = — f!g and 7 = Hq, -R(i) = exp(Oo*)- (b) Corresponding contour plot. The other paramete 
are a = 1, b = 0.5, c = 1.0, f^o = 0.1 in Eq. (43). 



|$-(z,t)| 



exp [i7oi] • 
(44) 



The exponential growth of the dark solitary wave due to the presence of gain is illustrated in Figs. 9(a) and 
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Fig. 9. (a) Dark soliton for ^■^(t) = — f^Q and 7 = Cqi = exp(f2ot). (b) corresponding contour plot. The other parameters 
are /3 = 1/^2, 6 = 0.5, c = 1.0, = 0.1 in Eq. (44). 

9(b), where the width decreases exponentially in time. 

5.2. Case 2: R = exp(-17ot) 

Next, in the case of i? = exp(— iloi) similar type of stable solitary waves and growth/decay of solitary 
waves can be identified for suitable forms of R{t) and 7(t). Also one may note here that the width of the wave 
will be increasing exponentially with time (i.e, proportional to exp(ir2oi)) in contrast to the case discussed 
in Sec. 5.1. Since the details are similar we do not present them here separately. 

5.3. Case 3: R = sech(fioi + S) 

Next, we analyze the solutions of Eq. (7) for the case R = sech(Oot + 6) and with different forms of gain 
term 7. 
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Fig. 10. Bright soliton : for Q^{t) = -fig and 7 = 0, R{t) = scch{not + S) case, (b) The eorresponding contour plot. The other 
parameters are a = 1, b = 0.5, c = 1.0, Qq = 0.1 in Eq. (45). 

5.3.1. 7 = 0, decay of bright/dark soliton and two-soliton bound state: 

First we choose 7 = 0. Consequently, R takes the form R = sech(i7o^ + as deduced from Eqs. (11) and 
(33). The bright sohton solution of (7) for 7 = can be written as 



|<f>+(z,i) 



-\/ scch(51o^ + S) 



X sech 



^/2a 



tanh(l^ot + 5) 



(45) 



Figs. 10 shows the decaying nature of the bright solitary wave 




Fig. 11. (a) Dark sohton : for n'^{t) = -Q^ and 7 = 0, fl(t) = sech{Cot + <5) case, (b) The corresponding contour plot. The 
other parameters arc fS = I/V2, b = 0.5, c = 1.0, Qo = 0.1 in Eq. (46). 



In the present case, another interesting structure similar to a bright two-soliton bound state [shown in 
Fig. 12] occurs when the propagation term is absent in the above soliton solution (i.e, velocity is zero) 
when the parameters b and c in Eq. (45) are chosen as & = 0, c = . One may note that a similar kind of 
two-soliton like bound state has been recently identified by Atre, Panigrahi and Agarwal using a different 
approach [40]. 

Next, the dark soliton solution of (7) for the present case can be written as 
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Fig. 12. (a) Bright soliton bound state: for n^(t) = — and 7 = 0, R(t) = scch(Qot + S) case, (b) Tiie corresponding contour 
plot. The other parameters are a = I ,b = 0, c = 0,Qq = 0.1 in Eq. (45). 



X + 4/3^ tanh^ 



NaB 



sech(f7ot + 5)] z — 



B ^ , ^'^o-B „\ tanh(fiot + 5) 



(46) 

Figs. 11 shows the decaying nature of the dark solitary wave solution for b ^ and c ^ 0, where the width 
also increases exponentially in time. Similar to the case of bright soliton solution, a dark two-soliton like 
bound state can arise when 6 = 0, c = in Eq. (46). Figs. 13 shows such a dark two-soliton like bound state 
for 5 = 0, c = in Eq. (46). 
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Fig. 13. (a) Darli soliton bound state: for Q'^{t) = — f2g and 7 = 0, R{t) = sech(f2ot -I- 5) case, (b) The corresponding contour 
plot. The other parameters are 13 = I/V2, b = 0, c = 0,^0 = 0.1 in Eq. (46). 



5.3.2. 7 = f2o tanh [Vlot + 5], bright/dark soliton: 

Next let us choose the form of the gain term as 7 = fio tanh [Hot 4 
scattering length becomes R = sech^(iloi + 5), as obtained from Eqs. 
bright soliton solution of (7) can be written as 



5]. Then the form of the atomic 
(11) and (33). Consequently, the 
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(47) 
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Fig. 14. (a) Bright soliton: for Q'^{t) = — Hq and 7 = QotaBh{Qot + 5), R{t) = sech(S7ot + 5) case.(b) The corresponding 
contour plot. The other parameters are a = 1 ,b = 0.5, c = 1.0, f^o = 0.1 in Eq. (47). 
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Fig. 15. (a) Dark soliton: for f2^(t) = — t^Q and 7 = ta.nh{flot + 5) , R{t) = sech(S7of + <5) case.(b) The corresponding contour 
plot. The other parameters are /3 = 1/^2, b = 0.5, c = 1.0,^0 = 0.1 in Eq. (48). 

Fig. 14(a) shows the corresponding bright solitary wave. One may note that the ampHtudc of the wave is 
constant (A+ = ■^^^) over all time due to the balance between the gain and atomic scattering length but the 

width of the wave increases with respect to time (proportional to ■\/cosh(ilot + d)) as shown in Fig. 14(b). 
Next, the corresponding dark soliton solution of (7) can be written as 
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(48) 



which is shown in Figs. 15. Here the amplitude, A 
time. 



, remains constant while the width increases with 



5.3.3. 7 = ait — a2 sin(Ki), collapse and revival of the bright/dark solitary waves: 

Finally, we consider the case with periodically varying gain term of the form 7 = ait — 02 sin(Ki). In this 
case, one may obtain the expressions for the atomic scattering length from Eqs. (11) and (33) as 



R = exp 



-ai- 



cos(Ki) 



sech(r2o^ + <^)- 



(49) 



19 



Here again we observe the periodic collapse and revival of the bright and dark solitary waves, simi- 
lar to the case of the confining potential as shown in Fig. 2. However, the amplitude is decaying with 

COs(Kt) 



time ^ sech(17ot + S) x exp 
X exp 



) while in Fig. 2 it increases with time y'sec(f2o^ + 1^) 



, cos{nt) 

ai- — — 



2k 



). However, we do not present the further details as the qualitative features are similar 
to the case 4.2. 

In the above we have identified various interesting soliton solutions that exist in the case of time- 
independent harmonic potentials. A summary of these results are tabulated in Table 1 for easy reference. 
It is also of wide interest at present to look at the nature of the solutions in the presence of time varying 
harmonic trap potentials as experiments in BEC with time varying harmonic trap potentials are also pos- 
sible. In the following section, we consider certain interesting soliton solutions that exist in a select set of 
time-dependent traps. 

Table 1 

Nature of soliton solutions for different forms of trap potential Q'^{t), modulated scattering length R{t) and gain/loss term7(t) 
in the time-independent harmonic traps 





m 


7(t) 


Nature of Solution 




(i) 


sec(f2ot) 





periodic peaks in amplitude 




(ii) 


r ai 2 cos(Kt)l 

sec(S2otj exp 1 — a2 

L 2 K J 


ait — a2 sin(Kt) 


periodic collapse & revival with amplitude growth 




(iii) 


sec^(f!ot) 


— f2o tan(f!()t) 


snake like structure with constant amplitude 




(i) 


exp(2rJot) 




amplitude constant / width decreases 




(ii) 


constant 


Ho 


amplitude grows/width decreases 




(iii) 


exp(-2r2ot) 


Qo 


amplitude constant / width increases 




(iv) 


constant 


-Oo 


amplitude decreases/ width increases 




(v) 


sech(f2oi) 





amplitude decreases/width increases 
or bound state occurs 




(vi) 


sech^(r2ot) 


tanh(f2ot) 


amplitude constant / width increases 




(vii) 


sech(r2oi) exp(f2oi) 


Ha 

2 


amplitude decreases/ width increases 




(viii) 


sech(not) exp [ "21*2 aa"""'"''] 


ait — a2 sin(ret) 


periodic collapse & revival with decaying amplitude 



6. Bright and dark solitons in the time-dependent harmonic traps 



In this section, we consider certain interesting time-dependent trap potentials for which Eq. (7) exhibits 
soliton solutions. We again identify a variety of soliton solutions by suitably choosing the gain and atomic 
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scattering length. 



6.1. Double modulated periodic potential 



First let us choose the double modulated periodic potential, ri^(t), of the form 

n'^it) = -Qo (— [1 - cos(2Ai)] - Acos(Ai) 



Substituting the above form of f2^(t) into Eq. (19), one may obtain a particular solution for R as [see 
Table C.l] 



R = exp 



-no 
X 



eos(Ai) 



(51) 



6.1.1. 7 = — r2osin(At), periodic collapse and revival of the two-soliton bound state 

In this case, the soliton solution can be identified by choosing a periodically varying gain or loss term of 
the form 7 = — f2osin(At). Then from Eqs. (11) and (33) one can get the atomic scattering length as 



R = exp 



-2not 

A 



cos(At) 



(52) 



In this case an interesting type of periodic bound state bright soliton solution of (7) can be obtained (for 

6 = 0, c = 0) as 



I , , Nub 

\<P+{z,t)\= exp 

a± 

X sech 



-fto cos(At) 



A 



^NoB 



■ exp 



— Oo cos(A<) 



(53) 
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Fig. 16. (a) Bright soliton bound state: for ^{t) = — Qo(f^o/2{l — cos(2At)) — Acos(At)) and 7 = — Slosin(At), 
R(t) = cxp[ cos(At)] case, (b) The corresponding contour plot. The parameters are a = 1, b = 0, c = 0, fig = X = 0.1 in 
Eq. (53). 



As before, one may observe that the amplitude of the soliton is 



exp 



— Oq cos(Ai) 
A 



and width is 



proportional to exp " — - indicating periodic collapse and revival of the wave function. In the above, 

we see an additional feature in the form of a periodic collapse and revival of the bright two-soliton bound 
state of the matter wave for b ~ and c = 0. The bound state indicates the absence of propagation term 
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Fig. 17. (a) Dark soliton: for Q^{t) = -r2o(Oo/2(l - cos(2At)) - Acos(At)) and 7 = -Qq sin(At), R{t) = exp[^^ cos(At)] case, 
(b) The corresponding contour plot. The parameters are /3 = l/\/2, b = 0, c = 0, Qq = X = 0.1 in Eq. (54). 

(for 6 = and c = 0). However collapse and revival phenomena occur due to the fact that amplitude 



exp 



— ^^o cos(Ai) 
A 



is periodically varying with time. This is illustrated in Fig. 16. 



The corresponding periodic bound state dark soliton solution of (7) is given as 



Nob 

2l3a_ 



■ exp 



-fio cos(At) 



A 



y. Ic^ + 4/32 tanh-^ 



Nag 



■ exp 



-f^o cos(At) 
A 



Fig. 17 shows the collapse and revival of the dark two-soliton bound state. 



6.2. Step-wise time- dependent trap potential 



(54) 



Next, let us consider a step-function modulation potential, Vt{t) of the form 



1 - tanh ( 



(55) 



Using the above form of VL{t) in the integrability condition (19), we get a particular solution for i?, see 
Table C.l, 



i? = 1 + tanh ( 



(56) 



6.2.1. Case 1: -f{t) = [l - tanh (^i)] , bright/dark soliton 
If we choose the time-dependent gain or loss term as 



lit) - 



1 — tanh ( —t 
2 



(57) 



then from Eqs. (11) and (33) the form of the atomic scattering length R turns out to be 

^_ 4exp(2riot) 
[1 + cxp(rJot)]^' 

Consequently, the bright soliton solution of (7) for the time-dependent potential (55) with the above loss 
term (57) can be written as 



(58) 
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=^ sech 



NOE 



1 + tanh 



X <^ 20oi + 4 log 



cosh 



~2~ 



— 2 tanh 



N^a% b N^al c 



(59) 



The bright sohtary wave solution, in this case, is shown in Figs. 18(a) and 18(b). One may note that the am- 
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Fig. 18. (a) Bright soliton: for (i) = 
corresponding contour plot in the t - 
time. The parameters are a = 1.0, b - 



[1 - tanh(^i)] , R{t) = l + tanh(^t) and 7 = 
X plane and (c) Snapshots of soliton profile in the z 
-- 0.2, c = 0.2, Qq = 0.1 in Eq. (59). 



tanh(^t)] case, (b) The 
X plane at different instances of 



plitude of the wave remains constant (^+ = ^ITT^ while the width gradually decreases (~ l/y^ [l + tanh (%^) 
with time. 

It is also of interest to look at the solution of the original GP equation (1) with the trap potential (55) 
with time-dependent loss as in (57). Fig. 18(c) shows snapshots of the bright soliton solution in the z — x 
plane at different instances of time. One may note that the set of soliton profiles which are shown Fig. 18(c) 
resembles the one that is obtained in certain recent experiments [4,5]. However, these experiments were 
carried out for the formation of bright matter wave solitons in ultracold ^Li gas, with confining/expulsive 
time independent harmonic trap potentials. Here we have identified a similar kind of soliton solution in the 
step- wise modulated harmonic potential, which is experimentally feasible to verify in the future. In Ref [57], 
Xue studied the bright soliton solution without gain term for this case. 

Next, in this case the dark soliton solution can be written as 



|<i>-(z,i)| 



4/32 tanh^ 



Nob 



1 + tanh 



not 



no 



2f7ot + 41og 



cosh 



not 



- 2 tanh 



not 



(60) 
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The dark soliton solution is shown in Fig. 19. In this case also the amplitude of the wave remains constant 
(i^)^ while the width gradually decreases + tanh (%^)]) with time. 




Fig. 19. (a) Dark soliton; for U^it) = -Hg [l - tanh{^t)] , R{t) = 1 +tanh(^t) and 7 = -^[1 - tanh(^t)] case, (b) The 
corresponding contour plot. The parameters are /3 = l/\/2, b = 0.2, c = 0.2, f!o = O.lin Eq. (60). 

6.2.2. Case ^; 7 = ^ > 0, growth of bright/dark soliton 

Finally, we consider a positive gain term of the form 7 = ^ > 0. For this choice of 7 Eqs. (11) and (33) 
yield the condition on the the atomic scattering length as i? = sech ( ) . In this case we see a growth of 
the bright /dark solitary wave similar to the case 5.1.2 with time independent potential. Here the amplitude 

^^^J sech [^] . exp [^] increases in time, while the width gradually decreases (-- + tanh (%^)]") 

with time. Since the details are similar to the case 5.1.2, we do not present them here separately. 



A summary of the results of time-dependent harmonic traps is tabulated in Table 2 for ready reference. 

Table 2 

Nature of soliton solutions for different forms of Q^{t), R{t) and ■y(t) in the time-dependent harmonic traps. 





R{t) 


7(t) 


Nature of Solution 


-Qo (^ [1 - cos(2At)] - Acos(At)) 


cxp cos(At)] 


— Qq sin Xt 


periodic collapse & revival of 
twro soliton like bound state 


-^[l-tanh(£^*)] 


4exp(2Qot) 


ai)^>o 


bright/dark soliton amplitude constant 

while width gradually decreases 
bright /dark soliton amplitude grows 


[1 -l-exp(noi)]^ 



7. Summary and conclusions 

In summary, we have investigated the exact bright and dark solitary wave solutions of an effective ID BEC 
by assuming that the interaction energy is much less than the kinetic energy in the transverse direction. In 
particular we have shown that the effective ID equation resulting from the GP equation can be transformed 
into the standard soliton (bright/dark) possessing, completely intcgrable ID NLS equation by effecting a 
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change of variables of the coordinates and the wave function. We have foeussed on both confining and 
expulsive harmonic trap potentials separately and treated the atomic scattering length, gain/loss term and 
trap frequency as experimental control parameters by modulating them as a function of time. Various 
types of bright and dark soliton solutions have been deduced in the context of effective ID Bose-Einstein 
condensation for different forms of unmodulated and modulated harmonic potentials by suitably tailoring 
the atomic scattering legnth and atom gain/loss term. First, we have observed the periodic oscillating 
solitons, collapse and revival of condensate and snake-like solitons for time independent confining harmonic 
potential. Next we have shown the existence of stable solitons, soliton growth and decay and formation of two- 
soliton bound state for time independent expulsive harmonic potential. However when the trap frequency 
is also modulated, we have shown the phenomena of collapse and revival of two-soliton like bound state 
formation of the condensate for double modulated periodic potential and bright and dark solitons for step- 
wise modulated potentials. Finally we have shown the snapshots of the bright soliton solution of 3D BEC in 
the z ~ X plane at different instances of time for step-wise modulated harmonic potential, which is similar 
to certain experimental [4, 5] results for time-independent potential. It appears that suitable experiments 
can be devised to identify these structures in BECs for time-dependent potentials. The experiments in BEC 
with these temporally modulated potentials are possible in future. Finally we will separately discuss the 
soliton interactions and the multi-soliton solutions of BECs with time varying parameters by extending our 
above analysis. 
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Appendix A. Validity for ID approximation for all cases in this paper 



The validity of the criteria for ID CP equation with time varying parameters as(i), T{t) and az{t) is given 
by Eq.(6). Now we substitute as{t) iii^, T{t) = a;^7(t), and a^X'>^) = a^y/\n{t)\ in Eq.(6) to get 



R{t)NaB^A^\ 
exp 



-f{t)dt 



« 1. (A.l) 



Then using the form R{t) = R{t) exp [/ 'y{t)dt], see Eq. (11), Eq. (A.l) for the validity of the ID approxi- 
mation can be rewritten as 

mNasVWm ^ (A 2) 

a± 

The above criterion should be satisfied for sufficiently long time for given form of ^{t) (which is related to 
az{t)) and R{t) (which is related as{t) and r(t)). The various forms of n(i) discussed in this paper are given 

in Tables 1 and 2. Fig. A.l shows the value of ^i^j^ ^^.^^ ^ foj. forms of R{t), n{t) studied 

in this paper. 

Fig. A. 1(a) shows the validity criteria of the ID CP equation for R{t) = sec{ilQt + 6) and ri^{t) = Oq. 
The validity criteria is satisfied for several milliseconds for R{t) = exp{^ot) and ^l'^{t) = —fig which is 
shown in Fig A. 1(b), which is sufficiently long for BEC experiments, See Sec. 5.1 for various forms of 
7(t). Figs. A. 1(c) and A. 1(d) show the validity criteria valid for all values of time by choosing R{t) = 
exp{—ilot) and R{t) = sech{V,ot + S) respectively for fl'^{t) = —^q, see Sec. 5.2 and 5.3. For Figs. (A.l 
a)- (A. Id), see also Table 1. Fig. A. 1(e) describes the validity criteria for R = exp [^cos(Ai)] and \n{t)\ = 

^/j — r^o [1 — cos(2Ai)] — Acos(At)) |, which shows the value jg periodically varying with 

time but within the validity regime. So the validity criteria is satisfied for all values of time. Fig. A. 1(f) 
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shows that the vaUdity criteria is valid for all values of time for R{t) = 1 + tanh (^t) and \^(t)\ = 
\J^- [l -tanh(^t)]. For Figs. (A.lc) and (A.lf), see also Table 2. 




5 10 15 20 25 30 5 10 15 20 25 30 

t (in ms) i (in ms) 




5 10 15 20 25 30 5 10 15 20 25 30 

t (in ms) / (in ms) 



Fig. A.l. Plot ( ._ (a) for R(t) = scc(not + <5) and = [see Sec. 4.1], (b) for R{t) = exp(Qot) and 

\n{t)\ = no [see Sec. 5.1], (c) for R{t) = cxp(-Qot) and \n(t)\ = Uq [see Sec. 5.2], (d) for R{t) = sech{not + S) and |Q(t)| = Ho 
[see Sec. 5.3], (e) for R = cxp cos(At)] and \n{t)\ = ^\ - f^o [1 - cos(2At)] - Acos{Ai)) | [see Sec. 6.1] and (f) for 

R(t) = 1 + tanh (^t) and |f2(t)| = \J \ ^ [l ~ tanh (^t)] | [sec Sec. 6.2]. The parameter values arc chosen as f!o = 0.1, 

5 = 0, A = 0.1, N = lO**, ax = 1.4/xm. For Figs, (a-d) sec also Table 1 and for (c) and (f) see Table 2 for more details. 



Appendix B. Solutions for the condition equations (16) for effective ID GP equation 
mapped onto the NLS equation 

Substituting the relation A(z, t) ~ r{t) exp[z6'(z, t)] from Eq. (17) into the condition Eqs. (16) and separate 
them into real and imaginary parts, we get the following 
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rt 



1 



^r ^e,,r ^ 0, 

^ -Z^ = 0, 



2 ^ 2 

Ft + e,F, = 0, 

F,, = 

Gt = ^ = i? 



From Eqs. (B.l) and (B.4), we obtain 

e 

F{z,t) =Fo{t)z + Fi{t), 



e{z,t) =—z^ + ai{t)z + a^{t), 
r 



(B.l 

(B.2 

(B.3 
(B.4 

(B.5 



(B.6 
(B.7 



where ai(t), a2{t) and Fi{t) are the functions of t to be determined. Substituting the above form of the 
Q{z,t) and F{z,t) into Eqs. (B.2) and (B.3) and cohecting the coefficient of the power z and then solving 
them, we get the following 



a2{t)=- 
Fo{t) 
Fi{t)=- 



al{t)dt 



d 


f2-l 




'2rt' 


di 


. r . 




r 



ai{t)Fo{t)dt 

2 

-n\t) = o 

where b and vq arc the integrating constants. From Eq. (B.5) gives 

G{t) =2 

Fo{t) ^rV2R 
Now comparing the Eqs. (B.IO) and (B.14), we get 

r = roV2^ 



r^Rdt 



(B.8 
(B.9 

(B.IO 

(B.ll 

(B.12 

(B.13 
(B.14 

(B.15 



Substituting the above form of r in to 0, ai, 0:2, Fq, Fi and G , we deduce the Eqs. (18). Then Eq. (19 
is obtained readily from the Eq. (B.12.) 



Appendix C. General solutions to Ricatti equation 



Here we describe some of the interesting solutions of the restrictive condition given by the Ricatti equa- 
tion (19). Substituting 



R = exp 



S{t)dt 



in Eq. (19), wc get a simplified form of the Ricatti equation as 



dS_ 

dt 



n\t) = 0. 



(C.l) 



The general solution of the above Ricatti equation (C.l) can be readily given as [71] 
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cxp [2 / Sodt] 



5 - / exp[2 / Sodt]dt\ ' 



(C.2) 



where Sq is a particular solution of equation (C.l) and c is an arbitrary constant. From the general solution 
of (C.l), one can write down the general solution of equation (19) using the relation R = exp[J S{t)dt]. 
There are numerous forms of J7^(t) for which explicit solutions for R(t) can be obtained [see for example 
Rcf. [71]]. Most important of them are given in the form of Table C.l 

Table C.l 

Explicit solutions of the Riccati equation (19) for different forms of ^l^{t). 





Form of n^{t) 


Physically interesting 


General solution of R 






particular solution of R 


( 


c: arbitrary constant) 


1. 


fig = constant 


sec{Qot) 




1 




cos(no*) 


— c sin(no i) 




2. 


— SIq = constant 


scch(not) 




1 








exp(±r2ot) 


cosh(nQ t) — c sinh(nQ£) 




3. 


-f[l-tanh 


1 + tanh {^t) 




-4c+(i+e"o')-4c(l+e''o*) iog(i+e''o') 


4. 


-f^O [1 - cos(2At)] + Acos{At)) 


exp ^ " cos{At) 


cxp 


cos{At) 




1-d J cxp]^ y° cos(Af)J dt 


5. 


r2oAexp(At) - nl exp(2At) 


exp (^exp(At)) 


exp(ii£!.cxp(At)) 


1 — c ExpIntcgralEi ^ cxp(At)^ 


6. 


3^0 - n^t^ 


1 f Slot^\ 


oxpl ^j— 


) 


_cxp(^^J 


c cxp(no 




7. 


flo - fi^t^ 


exp(lMi) 


cxp ^ 


' J 






l + y¥c Erfl(y^t) 




8. 






Complicated integral 


9. 


-2 


3t 


St 
















10. 




(2A+l)f^ 


Complicated integral 


(2A+l)+t2A+l ' 
where A = 


11. 


ClotX - f2o(rJo + A) coth2(At) 


[sinh(At)]""''^ 


Complicated integral 


12. 


noi2{l - 2coth2(not)) 


sinh(f!ot) 


sinh(f2Q i) 




l-2act + 


c sinh(2nQ £) 




13. 


3^0 tA - A2 - no(n(j + A)tanh2(At) 


[cosh{At)l"o/^ 


Complicated integral 


sinh(At) 
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Table C.l 
Continued. 



Form of n'^{t) 


Physically interesting R 
Particular solution R 


General solution of R 


14. 


2^2 sech2(rJot) 


coth(noi) 


cosh{nQ t) 
^-^ HQci) sinh(r2Q£) H-c cosh(£lQ£) 


15. 


SQoA - - Co(f^o + A) coth2(At) 


[sinh{At)l"o/^ 
cosh (At) 


Complicated integral 


16. 


-2^2 csech2(f1ot) 


tanh(noi) 


sinh(SlQi) 
(1 — Hoct) cosh (fiQ£) + csinh( Hot) 


17. 


Ho A - r2o (a + A) tanh^ (At) 


[cosh(At)]""/^ 


Complicated integral 


18. 


f7g(l - 2tanh2(f7ot)) 


cosh(not) 


cosh(S^Qi) 
l+2f2oc*+csinh(20ot) 


19. 


— S2q + S2oAsmn(Af) — S2q sinh (At) 


cxp 


|^£^sinh(At)) 


cxp(.j- sinh{A£)) 
1 — cj cxp(^ sinh(At)dt) 


20. 


-2Q2 [tanh2(r2oi) + coth2(r2o4)] 


cosh(not) sinh(noi) 


cosh{nQt 
l — 4Qo ct + 


sinh(noi) 
;sinh(4f2ot) 






21. 


-Qq + noAcos(At) + f^o cos^(At) 


cxp 


^^cos(At)) 


cxp 


f _ 


4^cos(At)^ 




1 — c J cxp 


y cos(At)dt 




22. 


-ng + f2oAsin(At) + sin'^{\t) 


exp 


(^sin(A<)) 


cxp 


t _ 


sin(At)' 




1 — c J cxp 


[^T^ sin(At)dt" 




23. 


QoA + r2o(A - Oo) tan2(At) 


sec(not)"o/^ 


Complicated integral 


24. 


A2 + sQoA + Ho (A - no) tan2(At) 


sin{Xt) 


Complicated integral 


25. 


A^ + 3^0 A + no (A - Ho) cot2(At) 


scc{nQt) 
sin(At)"o/^ 


Complicated integral 


26. 


no A + no (A - no) cot2(At) 


1 

sm{\t)"0/>- 


Complicated integral 


27. 


-2n2 [tan2(not) + cot2(not)] 


cos(not) sin(not) 


sin 

1 — 4^0 ct 


(200*) 

+ csin(4not) 







Appendix D. Physical meaning of symbols used in the text 

In order to make the various notations used in the text clear, we tabulate a list of symbols which we have 
used in the text with their physical meaning in Table D.l for ready reference. 
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